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Direct measurements of the spectra of extrasolar giant planets are the keys to determining 
their physical and chemical nature. The goal of theory is to provide the tools and context with 
which such data are understood. It is only by putting spectral observations through the sieve of 
theory that the promise of exoplanet research can be realized. With the new Spitzer and HST 
data of transiting "hot Jupiters," we have now dramatically entered the era of remote sensing. 
We are probing their atmospheric compositions and temperature profiles, are constraining their 
atmospheric dynamics, and are investigating their phase light curves. Soon, many non-transiting 
exoplanets with wide separations (analogs of Jupiter) will be imaged and their light curves and 
spectra measured. In this paper, we present the basic physics, chemistry, and spectroscopy 
necessary to model the current direct detections and to develop the more sophisticated theories 
for both close-in and wide-separation extrasolar giant planets that will be needed in the years to 
come as exoplanet research accelerates into its future. 



1. INTRODUCTION 

Our understanding of gas giant planets was informed 
for many decades by remote telescopic observations and in 
situ measurements of Jupiter and Saturn. These detailed 
investigations provided a fine-grained view of their atmo- 
spheric compositions, temperatures, dynamics, and cloud 
structures. However, they left us with a parochial view of 
the range of possible orbits, masses, and compositions that 
has now been shattered by the discovery of extrasolar giant 
planets (EGPs) in the hundreds. We have found gas giants 
in orbits from ^0.02 AU to many AU, with masses from 
below Neptune's to ~10 Mj , and around stars from M to 
F dwarfs. The corresponding stellar irradiation fluxes at the 
planet vary by a factor of ~10 5 , and this variation trans- 
lates into variations in atmosphere temperatures from ^100 
K to ^2500 K. With such a range of temperatures and of 
orbital distances, masses, and ages, atmospheres can have 
starkly different compositions, can be clear or cloudy, and 
can evince dramatic day-night contrasts. 

One must distinguish imaging of the planet itself by sep- 
arating the light of planet and star, something that can cur- 
rently be contemplated only for wide-separation planets, 
from measurements of the summed light when the orbit is 
tight and the planet can not be separately imaged. In the 
latter case, the planet's light can be a non-trivial fraction 
of the total, particularly in the infrared. When transiting, 
such hot Jupiter systems provide an unprecedented oppor- 
tunty to measure the planet's emissions by the difference in 
the summed light of planet and star in and out of secondary 
eclipse and by the phase variation of that sum. Generally, 
the star itself will not vary with the period of the planet's 
orbit. Moreover, in a complementary, but different, fash- 



ion, the wavelength dependence of the transit depth is now 
being used to probe the composition of the planet's atmo- 
sphere near the terminators. The extrasolar giant planets, 
by dint of their mass and luminosity, have been the first dis- 
covered, and will serve as stepping stones to the extrasolar 
terrestrial planets. 

To understand in physical detail the growing bestiary 
of EGPs requires chemistry to determine compositions, 
molecular and atomic spectroscopy to derive opacities, ra- 
diative transfer to predict spectra, hydrodynamics to con- 
strain atmospheric dynamics and heat redistribution, and 
cloud physics. In short, global 3D radiation-hydrodynamic 
general circulation models (GCMs) with multi-spectral, 
multi-angle, and non-equilibrium chemistry and kinetics 
will be needed. We are not there yet, but basic treatments 
have emerged that allow us to interpret and constrain day- 
night differences, profiles, molecular compositions, and 
phase light curves. 

In this chapter we lay out some of the basic elements of 
any theoretical treatment of the atmospheres, spectra, and 
light curves of EGPs. This theory provides the necessary 
underpinnings for any progress in EGP studies, a subject 
that is engaging an increasing fraction of the world's astro- 
nomical and planetary science communities. In §2.1, we 
summarize the techniques for calculating molecular abun- 
dances. We follow in §2.2 with an explication of general 
methods for assembling opacity tables. Section 2.3 touches 
on Rayleigh scattering, and then we continue in §2.4 with a 
tutorial on albedos and phase functions. In §2.5, we explain 
the nature of the transit radius. Section 2.7 contains a very 
useful analytic model for the atmospheric thermal profile of 
EGPs, which is a generalization for irradiated atmospheres 
of the classic Milne problem. This model incorporates a 



1 



condition for thermal inversions. Then, in §3 we summarize 
lessons learned and knowledge gained from the decades- 
long study of Jupiter and Saturn. This includes discussions 
of their spectra, cloud layers, temperatures, and composi- 
tions. Having set the stage, we review in §4.1 the general 
chemistry and atmospheric character of EGPs as a func- 
tion of orbital distance and age. This subsection includes 
a diversion into the putative evolution of Jupiter itself. We 
follow this in §4.2 with a few paragraphs on theoretical 
EGP planet/star flux ratios as a function of wavelength and 
distance, focusing on wide-separation (> 0.2 AU) EGPs. 
Then, in §4.3 we present highlights from recent campaigns 
of direct detection of transiting EGP atmospheres, with an 
emphasis on secondary eclipse measurements, the compo- 
sitions inferred, atmospheric temperatures, and thermal in- 
versions. Finally, in §5 we list some of the outstanding open 
issues and future prospects in the study of EGP atmospheres 
and spectra. 

2. MODELING CONCEPTS AND EQUATIONS 

In this section, we present the core ingredients neces- 
sary to construct theories of the atmospheres and spectra 
of EGPs. These includes general chemistry, opacities, and 
simple formalisms for the calculation of albedos and phase 
light curves. We discuss the concept of "transit radius" and 
include an analytic theory for the temperature profiles of ir- 
radiated exoplanets. The resulting approximate equations 
make clear the key role of opacity and its wavelength de- 
pendence in determining the character of EGP thermal pro- 
files, in particular in creating thermal inversions when they 
arise. Such inversions have been inferred for many of the 
hot Jupiters seen in secondary eclipse and are emerging as 
one of the most exciting and puzzling features of current 
exoplanet research. 

2.1 Calculation of Atomic and Molecular Abundances 
in Chemical Equilibrium 

Before any opacities or atmopsheric models can be cal- 
culated, the abundances of a mixture of a large number of 
species have to be determined for the given temperature and 
pressure. The assumption of chemical equilibrium is a good 
starting point, though non-equilibrium kinetics may play a 
role. Despite this, we present here a straightforward dis- 
cussion of such calculations. Much of this presentation on 
abundances and that of §2.2 on opacities can be found in 
Burrows & Sharp (1999) or Sharp & Burrows (2007), to 
which the reader is referred. In addition, there are excellent 
reviews and/or accounts in Lodders (1999) Fegley & Lod- 
ders (1994,2001), Lodders & Fegley (2002), and Sharp & 
Huebner(1990). 

For a given temperature, pressure, and composition, the 
equilibrium abundances of the various species can be deter- 
mined by minimizing the total Gibbs free energy of the sys- 
tem. This requires a knowledge of the free energy of each 
species as a function of temperature, which is normally ob- 
tained from thermodynamic data. At the temperatures for 



which data are tabulated, least-square fits can be made for a 
set of polynomials whose highest order is given by 



AG pl (T) = aT- 1 + b + cT + dT 2 + eT 3 



(1) 



where a, b, c, d, and e are fitted coefficients, and AG P i(T) 
is the fitted Gibbs free energy of formation at temperature 
T of species i in phase p. The polynomials are evaluated 
at the tabulated points and the deviations from the tabulated 
values are obtained. 

In performing the calculation for a particular tempera- 
ture, pressure, and composition, the Gibbs free energy for 
each species is obtained from the database using the fitted 
coefficients at the temperature required, then the total free 
energy of the system is minimized to obtain the abundances 
of the gas-phase species, together with any condensates. 
The total dimensionless free energy is given by 
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where R is the gas constant (kgN^), and for the first sum 
for the gas phase with p = 1 , P is the total pressure in 
atmospheres, N is the number of moles, m is the number 
of species, nu is the number of moles of species i, and 
AGh(T) is the corresponding free energy of that species. 
The second sum is over the s condensed phases, which may 
include multiple phases of the same species, but except at a 
phase boundary, only one phase of a particular species in a 
condensed form is assumed present at any time, though one 
may consider solid or liquid solutions. Consequently, n p i is 
the number of moles of a condensed species and AG p \{T) 
is the free energy of that species. Since there is only one 
species per phase, for convenience we generally set i equal 
to 1. 

A subset of 30 gas-phase species out of nearly 350 gas- 
phase species are usually the most important and selected 
for detailed treatment. The species are the neutral atoms: 
H, He, Li, Na, K, Rb, Cs, Al, Ca, and Fe, the ions: e", 
H+, and H", the metal hydrides: MgH, CaH, FeH, CrH 
and TiH, with the remaining molecules being H 2 , N 2 , CO, 
SiO, TiO, VO, CaOH, H 2 0, H 2 S, NH 3 , PH 3 , and CH 4 . 

Figure Q] depicts a representative result (here at 1 atmo- 
sphere pressure and solar elemental abundances) for the 
temperature dependence of the equilibrium mixing ratios 
of water, methane, carbon monoxide, hydrogen sulfide, 
phosphene, molecular nitrogen, ammonia, TiO, and VO. 
The last two species may or may not play a role in EGP 
atmospheres, whereas the others certainly do. 

2.2 Calculation of Atomic and Molecular Opacities 

The calculation of the absorption cross sections and 
opacities of molecules is made more difficult than that for 



2 



atoms by the substantially larger number of transitions and 
levels involved. Polyatomic species can have hundreds of 
millions, even billions, of vibrational and rotational lines, 
multiple electronic states, and a complicating mix of iso- 
topes. Since it is not possible to measure with precision 
many transitions to determine their oscillator strengths or 
Einstein A coefficients, ab initio calculations using quan- 
tum chemical techniques are frequently necessary. Such 
calculations can be, and frequently are, calibrated with only 
a few measurements at selected wavelengths, but the exper- 
imental determination of the quantum numbers of the upper 
and lower states of even a given measured transition can 
be ambiguous. Moreover, particularly for hot Jupiters, the 
high temperatures experienced require a knowledge of ab- 
sorption transitions from excited states, the so-called "hot 
bands," for which there are rarely measurements. For in- 
stance, the methane hot bands and some of the hot bands of 
water are completely unconstrained by current experiment, 
despite the fact that both water and methane are important 
greenhouse gases in the Earth's atmosphere. 

The situation is made even more difficult by the almost 
complete absence of calculations or measurements of line 
broadening coefficients or line profiles. The upshot is that 
theorists rely on imperfect and by-and-large uncalibrated 
compilations of theoretical line lists, and very approximate 
theories for line broadening. Some molecules are done bet- 
ter than others, but even for atoms, for which theory and 
measurement are rather better, line shapes are not well con- 
strained. In the context of EGPs, this is particularly relevant 
for the alkali metals. 

Despite these drawbacks, a rather sophisticated and ex- 
tensive database of opacities for the constituents of EGPs 
has been assembled. In this subsection, taken in part from 
the review paper by Sharp & Burrows (2007), we summa- 
rize the techniques and methodolgies needed to derive these 
opacities and use cgs units for specificity. Another useful 
review is that of Freedman, Marley, & Lodders (2008). 

The calculation of the line strengths for each line of each 
species depends on the data available for the species being 
considered, so different methods have to be used. In order 
to reduce the chances of errors with input data in different 
forms, it is recommended to convert, if necessary, all the 
line strengths into the same uniform system, with the best 
being integrated line strengths in cm 2 s _1 species -1 . These 
depend only on the temperature. The lines should then be 
broadened into a profile which is dependent on the pres- 
sure, then the absorption in cm 2 species _1 across the profile 
should be computed, summing the contributions from any 
overlapping profiles. The absorption for each species ob- 
tained in this manner depends on the temperature and pres- 
sure. The total opacity of the gas is obtained by summing 
the individual contributions weighted by the corresponding 
number densities (in cm -3 ) for each of the species, yielding 
the total volume opacity in cm 2 cm~ 3 , i.e. cm -1 . However, 
the total mass opacity in cm 2 g _1 is usually the required re- 
sult, and is obtained by dividing the volume opacity by the 
gas mass density. 



In its most general LTE (Local Thermodynamic Equilib- 
rium) form, the integrated strength S of a spectral line in 
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where gi and /y are the statistical weight of the i energy 
level and the oscillator strength for a transition from that 
level to a higher level j, Fi, and Fj are the term values (ex- 
citation energies) in cm -1 of the i th and j th levels partici- 
pating in the transition, and Q(T) is the partition function of 
the species at some temperature T. The other symbols have 
their usual meanings. Note that the first term in eq. (01 gives 
the line strength in cm 2 s _1 absorber -1 , the next term with 
the Boltzmann factor and the partition function converts this 
to the required line strength, and the last term is the stimu- 
lated emission correction factor, where Fj — Fi is the tran- 
sition frequency in wavenumbers, i.e. v in cm -1 . Although 
monochromatic opacities are frequently displayed as func- 
tions of wavelength, it is recommended that all opacity cal- 
culations be performed internally in wavenumbers, even if 
some of the input data are given in wavelengths, since most 
molecular spectroscopic constants and energy levels are ex- 
pressed in cm -1 , and adopting a uniform system of units re- 
duces the chances of error. Note that some of the data avail- 
able are not expressed in the form of oscillator strengths 
and statistical weights, collectively given as (//-values, but 
in other forms that must be converted to the required lines 
strengths. 

The general expression for calculating the partition func- 
tion is given by 



Q(T) 



(4) 



where the summation is performed over the first n levels, 
whose contributions are required at the highest tempera- 
tures of interest. The term value of the lowest level F\, 
i.e. the gound state, is zero by construction. 

Figure |2]portrays a representative set of absorption cross 
sections as a function of wavelength in the infrared for the 
vibration-rotation transitions of H2O, NH3, andCITt. Many 
of the most important absorption bands in EGP atmospheres 
are shown and these were calculated using the formalism 
described in this subsection. 

2.3 Rayleigh Scattering 

Rayleigh scattering is a conservative scattering process 
by atoms and molecules. Although strong in the ultravi- 
olet/blue, the scattering cross sections quickly fade toward 
the red region of the spectrum (oc A~ 4 ). Rayleigh scattering 
has little effect on the spectra of isolated brown dwarfs, but 
irradiated EGPs reflect a non-zero fraction of the incident 
intensity. 
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The Rayleigh scattering cross sections are derived from 
polarizabilities, which are in turn derived from refractive 
indices. The refractive indices are readily available at 5893 
A (Na D) and are assumed not to vary strongly with wave- 
length. The Rayleigh cross sections are derived via, 
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where n is the index of refraction, k is the wavenum- 
ber (2tt/X), and Lq is Loschmidt's number, the number 
of molecules per cubic centimeter at STP (= 2.687 x 10 19 ). 
Given the strong inverse dependence on wavelength of eq. 
(0, Rayleigh scattering is most pronounced in the blue and 
ultraviolet and is ultimately responsible for the Earth's blue 
sky. 

2.4 Albedos and Phase Functions 

Objects in the solar system, such as planets, asteroids, 
and moons, are seen and studied in reflected solar light. 
The brightness of the reflection depends upon the orbital 
distance, the stellar flux, the reflectivity of the object, the 
detector angle (the "phase angle"), and the object's radius. 
The reflectivity, in the guise of an "albedo" (defined below), 
bears the stamp of the composition of its surface and/or at- 
mosphere, and its wavelength dependence is a distinctive 
and discriminating signature. Solar-system objects are too 
cold to emit much in the optical, where the Sun is brightest, 
but emit in the mid-infrared, where the peak of a ~40-800 
K black body resides. Hence, there is a simple and obvi- 
ous separation in their spectra between reflection and emis- 
sion bumps that allows an unambiguous definition of the 
albedo and its interpretation as a dimensionless reflectivity 
bounded by a value of "one" (however, see below), the latter 
implying full reflection and no absorption. 

However, some EGPs, the close-in and transiting vari- 
ants, are so near their central stars that their surface and 
atmospheric temperatures can be quite large (^1000-2500 
K). Hot planets (due to either proximity or youth) can be 
self-luminous in the near infrared (and even in the optical). 
As a result, the reflection and emission components can 
overlap in wavelength space, and are not so cleanly sepa- 
rated, as they are for solar-system objects. The upshot is 
that the albedo and "reflectivity" might be misnomers, par- 
ticularly in the near-IR. Nevertheless, the planet/star flux 
ratio as a function of wavelength is an important probe of 
EGP atmospheres and has to date been used with profit to 
diagnose their thermal and compositional character. This is 
true despite the fact that the associated albedos could be far 
above one at some wavelength^- With this caveat in mind 
and the traditional interpretation intact for planets at greater 
orbital distances and older ages similar to the Jupiter/Sun 
pair, we proceed to develop the formalism by which the 



1 Note that if the complete radiative transfer solution with stellar irradiation 
is derived, the concept of an albedo is redundant and unnecessary. 



planet/star flux ratio is calculated, the albedo is defined, and 
why they are important. In the process we distinguish the 
geometric albedo (A g ), the spherical albedo (A s ), and the 
Bond albedo (As), and connect them to the planet/star flux 
ratio, Fp/F*. We also introduce the "phase function," by 
which the observed light curve is described. Much of the 
development below is taken from papers by Sudarsky, Bur- 
rows, & Pinto (2000), Sudarsky, Burrows, Hubeny, & Li 
(2005), and Burrows, Sudarsky, & Hubeny (2004). The pa- 
pers by Marley et al. (1999) and Burrows, Ibgui, & Hubeny 
(2008) are also good resources on this general topic. 

Planetary phase is a function of the observer-planet-star 
orientation, and the angle whose vertex lies at the planet 
is known as the phase angle (a). The formalism for the 
computation of planetary brightness as a function of phase 
angle has been presented by numerous authors. Following 
Sobolev (1975), one can relate the planetary latitude (ip) 
and longitude (£) to the cosine of the angle of incident radi- 
ation (/.to) and the cosine of the angle of emergent radiation 
(p) at each point on the planet's surface: 



and 



p = cos V> cos(a — £) 



p = cos ip cos £ . 



(6) 
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where latitude is measured from the orbital plane and lon- 
gitude is measured from the observer's line of sight. The 
phase angle is then, 



a = cos 
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where 4> is the azimuthal angle between the incident and 
emergent radiation at a point on the planet's surface. The 
emergent intensity from a given planetary latitude and lon- 
gitude is 

I(p,p ,<j>) = p Sp(p,p ,4>), (9) 

where the incident flux on a patch of the planet's surface is 
npoS, and p(p, po, (j>) is the reflection coefficient. In order 
to compute the energy reflected off the entire planet, one 
must integrate over the surface of the planet. For a given 
planetary phase, the energy per time per unit area per unit 
solid angle received by an observer is 

R 2 r /2 

E(a) = 2S-£ / cosfa - cos(£)d£ x (10) 
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where R p is the planet's radius and d is the distance to 
the observer. This quantity is related to the geometric 
albedo (A g ), the reflectivity of an object at full phase 
(a = 0) relative to that of a perfect Lambert disk (for which 
p(p, po,4>) = 1) of the same radius under the same incident 
flux, by 

E(0)d 2 
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where E(0) is E(a) at a = 0. A planet in orbit about its 
central star displays a range of phases, and the planet/star 
flux ratio is given by 



*(o), 



(12) 



where $(a) is the classical phase function, which is equal 
to E(a)/E(Q), R p is the planet's radius, and a is its orbital 
distance. This formula is one of the core relationships in the 
study of irradiated and reflecting EGPs. 

<fr(a) is normalized to be 1.0 at full face, thereby defin- 
ing the geometric albedo, and is a decreasing function of 
a. For Lambert reflection, an incident ray on a planetary 
patch emerges uniformly over the exit hemisphere, A g is 
2/3 for purely scattering atmospheres, and $(a) is given by 
the formula: 



$(a) 



sin(a) + (tt — a) cos(a) 



(13) 



However, EGP atmospheres are absorbing and the anisotropy 
of the single scattering phase function for grains, droplets, 
or molecules results in non-Lambertian behavior. For in- 
stance, back-scattering off cloud particles can introduce an 
"opposition" effect for which the planet appears anoma- 
lously bright at small a. Figure [3]provides some theoretical 
EGP phase functions taken from Sudarsky et al. (2005) 
in which this effect is clearly seen. Figure |4] depicts the 
corresponding phase curves for some the solar system ob- 
jects. These two figures together suggest a likely range for 
exoplanets. 

Both A g and $(a) are functions of wavelength, but the 
wavelength-dependence of A g is the most extreme. In fact, 
for cloud-free atmospheres, due to strong absorption by 
molecular bands, A g can be as low as 0.03, making such 
objects very "black." Rayleigh scattering serves to support 
A g , but mostly in the blue and UV, where various exotic 
trace molecules can decrease it. The presence of clouds in- 
creases A g significantly. For instance, at 0.48 fim , Jupiter's 
geometric albedo is —0.46 and Saturn's is 0.39 (Karkoschka 
1999). However, for orbital distances less than 1.5 AU, 
we expect the atmospheres of most EGPs to be clear. The 
albedo is correspondingly low. As a consequence, the theo- 
retical albedo is very non-monotonic with distance, ranging 
in the visible from perhaps —0.3 at 0.05 AU, to —0.05 at 
0.2 AU, to -0.4 at 4 AU, to -0.7 at 15 AU. In the vis- 
ible (—0.55 /im), the geometric albedo for a hot Jupiter 
is severely suppressed by Na-D at 0.589 /im . Due to a 
methane feature, the geometric albedo can vary from 0.05 
at —0.6 /mi to —0.4 at 0.625 /im. Hence, variations with 
wavelength and with orbital distance by factors of 2 to 10 
are expected. 

$(a) and A g must be calculated or measured, but the 
sole dependence of <£>(a) on a belies the complications in- 
troduced by an orbit's inclination angle (i), eccentricity (e), 
and argument of periastron (co). Along with the period (P) 
and an arbitrary zero of time, these are most of the so-called 



Keplerian elements of an orbit. Figures 5 and 6 in Chapter 2 
of this volume define these orientational and orbital param- 
eters. In the plane of the orbit, the angle between the planet 
and the periastron/periapse (point of closest approach to the 
star) at the star is 9. In celestial mechanics, 9 is the so-called 
"true anomaly." For an edge-on orbit (i = 90°), and one for 
which the line of nodes is perpendicular to the line of sight 
(longitude of the ascending node, 17, equals 90°) and paral- 
lel to the star-periapse line (u> = 0°),6 is complementary to 
a (a = 90° - 9). As a result, 9 = 0° at a = 90° (great- 
est elongation) and increases with time. Also, for such an 
edge-on orbit, a = 0° at superior conjunction. In general, 

cos(a) = s'm(9 + u) sin(i) sin(17) — cos(fi) cos(# + u>) . 

(14) 

Since it is common to define the observing coordinate sys- 
tem such that 17 = 90°, we have the simpler formula: 



cos(a) = sin(# + oj) sin(i). 



(15) 



In order to produce a model light curve for a planet or- 
biting its central star, one must relate the planet's orbital 
angle (6, the true anomaly), as measured from periapse, to 
the time (t) in the planet's orbit: 



t{9) = 



.(l_ e 2)l/2p 
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2(1 -e 2 )- 1 ' 2 tan" 



(l-e 2 ) 1/2 tan(0/2) 



where P is the orbital period and e is the eccentricity. By 
combining eq. ( fTSl and ( fT5l ), we derive the exact phase of 
any orbit at any time. 

For a circular orbit, R is equal to the semi-major axis 
(a). However, a planet in an eccentric orbit can experience 
significant variation in R, and, therefore, stellar irradiation 
(by a factor of (y^f) 2 ). For example, if e = 0.3, the stellar 
flux varies by —3.5 along its orbit. For e = 0.6, this varia- 
tion is a factor of 16! Such eccentricities are by no means 
rare in the sample of known EGPs. Therefore, in response 
to a changing stellar flux it is possible for the composition 
of an EGP atmosphere to change significantly during its or- 
bit, for clouds to appear and disappear, and for there to be 
lags in the accommodation of a planet's atmosphere to a 
varying irradiation regime. Ignoring the latter, eqs. ( TT2T > 
and ( TBl i can be combined with $(a) and the standard Ke- 
plerian formula connecting 9 and time for an orbit with a 
given P and e to derive an EGP's light curve as a function 
of wavelength, i, e, 17, u, and time. Depending upon orien- 
tation and eccentricity, the brightness of an EGP can vary 
in its orbit not at all (for a face-on EGP in a circular or- 
bit) or quite dramatically (e.g., for highly eccentric orbits at 
high inclination angles). Since astrometric measurements 
of stellar wobble induced by EGPs can yield the entire orbit 
(including inclination), data from the Space Interferometry 
Mission (SIM) (Unwin & Shao 2000) or Gaia (Perryman 
2003) could provide important supplementary data to aid in 
the interpretation of direct detections of EGPs. 
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The spherical albedo is the fraction of incident light re- 
flected by a sphere at all angles. For a theoretical object 
with no absorptive opacity, all incident radiation is scat- 
tered, resulting in a spherical albedo of unity. The spher- 
ical albedo is related to the geometric albedo by A s = qA g , 
where q is the phase integral: 



(17) 



q = 2 $(a) sin ada . 



For isotropic surface reflection (Lambert reflection) q = |, 
while for pure Rayleigh scattering q = | . Although not 
written explicitly, all of the above quantities are functions 
of wavelength. 

Van de Hulst (1974) derived a solution for the spheri- 
cal albedo of a planet covered with a semi-infinite homo- 
geneous cloud layer. Given a single-scattering albedo of 
a (= (J scat I '(Ttotai) and a scattering asymmetry factor of 
g =< cos 9 > (the average cosine of the scattering angle), 
van de Hulst's expression for the spherical albedo of such 
an atmosphere is 



A s 



(1 -0.139s)(l-s) 
1 + 1.170s 



where 



1 



1/2 



(18) 



(19) 



1-0-0 

The Henyey-Greenstein single-scattering phase func 



tion, 



p(0) = 



9 



(20) 



(l + .g 2 -2 5 cos<9) 3 / 2 ' 

is frequently used as a fit to the overall scattering phase 
function, where again g =< cos 9 >. Other phase functions 
can be used, but their specific angular dependence has been 
found to be less important than the value of the integral, 
g, itself. In addition, Rayleigh and cloud particle scatter- 
ing are both likely to result in significant polarization of the 
reflected light from an EGP (Seager, Whitney, & Sasselov 
2000). The degree of polarization will depend strongly on 
orbital phase angle and wavelength, and can reach many 
tens of percent, However, polarization may be difficult to 
measure. To date, there is no credible evidence for polar- 
ized light from any EGP. The degree of polarization in the 
optical and UV is expected to be largest. 

The important Bond albedo, Ab, is the ratio of the total 
reflected and total incident powers. It is obtained by weight- 
ing the spherical albedo by the spectrum of the illuminating 
source and integrating over all wavelengths: 
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where the A subscript signifies that the incident intensity 
varies with wavelength. 

With the Bond albedo, one can make a crude estimate of 
the "effective temperature" of the planet's emission compo- 
nent in response to irradiation. Under the assumption that 



the incident total power is equal to the emitted power, one 
derives: 
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where T* is the stellar effective temperature, a is the orbital 
distance, i?* is the stellar radius, and / is a measure of the 
degree of heat redistribution around the planet. It is equal 
to 1/4 if the reradiation is isotropic, and 1/2 if the planet 
reemits only on the day side and uniformly. Note that T e s 
is independent of the planet's radius. 

Eq. d22l is useful in estimating the temperatures 
achieved by radiating planets, but has conceptual limita- 
tions. First, it presumes that the planet is not self-luminous 
and radiating the residual heat of formation in its core. This 
assumption is not true for young and/or massive EGPs. It 
is not true of Jupiter. Second, it encourages the notion 
that the reflection and emission peaks are well-separated. 
This is not correct for the hot Jupiters, for which the two 
components overlap and merge. Third, in radiative trans- 
fer theory and equilibrium, the true flux from a planet is 
the net flux. Without internal heat sources, this is zero for 
irradiated planets. Finally, eq. d22l is often used to de- 
termine the temperature of an atmosphere. However, EGP 
atmospheres have temperature profiles. In their radiative 
zones the temperatures can vary by factors of ~3, as can 
the effective photospheric temperatures in the near- and 
mid-IR. Hence, one should employ eq. d22"b to obtain atmo- 
spheric temperatures only when very approximate numbers 
are desired. Importantly, calculating the irradiated planet's 
spectrum with a T e g derived using eq. d22l) results in large 
errors across the entire wavelength range that can severely 
compromise the interpretation of data. 

2.5 The Transit Radius 

A transiting planet reveals its radius (i? p ) by the magni- 
tude of the diminution in the stellar light during the planet's 
traverse of the stellar disk. This is the primary eclipse. In 
fact, it is the ratio of the planet and star radii that is most 
directly measured, so a knowledge of the star's radius is 
central to extracting this important quantity. With a radius 
and a mass (for transiting planets, the inclination must be 
near 90° and the inclination degeneracy is broken), we can 
compare with theories of the planet's physical structure and 
evolution. However, since the effective edge of the plan- 
etary disk is determined by the opacity of the atmosphere 
at the wavelength of observation, the radius of a gas giant 
is wavelength-dependent. Importantly, the variation with 
wavelength of the measured radius can serve as an ersatz 
atmospheric "spectrum." From this spectrum, one can de- 
termine the atmosphere's constituent atoms and molecules. 
Specifically, the apparent radius is larger at wavelengths 
for which the opacity is larger and smaller at wavelengths 
for which the opacity is smaller. For example, the planet 
should be larger in the sodium D line at ~5890 A(as in fact 
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was found for HD 209245b) than just outside it. It should 
be larger near the absorption peaks of the water spectrum 
for those planets with atmospheric water than in the corre- 
sponding troughs. 

However, the "transit radius" is not the same as the clas- 
sical photospheric radius of an atmosphere. The latter is 
determined by the depth in the atmosphere where t v = 2/3 
in the radial direction. The transit radius is where this same 
condition obtains along the chord from the star, perpendicu- 
lar to the radius to the center. At the spherical radius where 
a light beam experiences r v — 2/3, r„ along the chord 
can be much larger. Therefore, to achieve the t v — 2/3 
condition along the chord pushes the transit radius (also 
referred to as an "impact parameter") to larger values. It 
is only at this greater altitude and lower pressure that the 
chord optical depth is <~2/3. An important difference be- 
tween the transit radius spectrum and the actual "emission" 
spectrum from a classical photosphere (such as is relevant 
at secondary eclipse) is that in the latter case, if the atmo- 
sphere were isothermal the spectrum would be a black body 
and there would be no composition information. However, 
the transit radius spectrum always manifests the wavelength 
dependence of the opacities of its atmospheric constituents, 
even if the atmosphere were isothermal. This difference 
can be exploited to maximize the scientific return from the 
study of a given transiting EGP. In sum, if the atmosphere is 
extended and if the monochromatic opacity at the measure- 
ment wavelength is large, the transit radius and the pho- 
tosphere radius can be rather different and the distinction 
should always be kept in mind. 

A difficulty in interpreting transit radius spectra is that 
one is probing the planetary limb, the terminator. This 
means that when comparing with data models must incor- 
porate profiles on both the day side and night side, and at 
both the equator and poles. In particular the day side and 
the night side can be at different temperatures and have dif- 
ferent compositions. This complication is not always ap- 
preciated. 

One can estimate the magnitude of the excess of the 
transit radius over the photospheric radius using a sim- 
ple exponential atmosphere (see also Fortney 2005). The 
wavelength-dependent optical depth, r c h or d, along a chord 
followed by the stellar beam through the planet's upper at- 
mosphere, is approximately: 



|, this yields 
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This excess can be from <~1% to ^10%, depending upon 
the wavelength, temperature, gravity, and deviations from 
a strictly exponential profile. It is smallest for high-gravity 
planets at larger orbital distances, whereas in the UV near 
Lyman-a, and with a planetary wind, this excess for HD 
209458b is measured to be a factor of 2—3. 

At times, the transit radius spectrum is called the "trans- 
mission spectrum." Since what is measured is the effective 
area of the planet, not the spectrum of light transmitted 
through the planetary limb, this is a slight misnomer. One is 
measuring the atmospheric edge position (relative to some 
zero-point) as a function of wavelength, and not the light 
transmitted through the finite extent of the atmosphere. To 
properly do the latter would require a resolved image of the 
planet that distinguishes the atmosphere from the opaque 
central planetary disk. Note also that when the atmospheric 
opacity is low (and, hence, light is more easily "transmit- 
ted"), the dip in the stellar flux (what is actually measured) 
is smaller, not larger. Nevertheless, this rather pedantic 
point does not inhibit one from profiting from "transmis- 
sion spectrum" measurements. 

2.6 Analytic Model for the Temperature Profile of an 
Irradiated Planet 

One can derive a model for the temperature profile of an 
irradiated EGP atmosphere that is a generalization of the 
classical Milne problem for an isolated atmosphere. This 
model can incorporate the difference between the opacity 
to the insolating ("optical") and emitting ("infrared") radi- 
ation streams, contains a theory for stratospheres and ther- 
mal inversions, and a condition for their emergence. The 
mathematical development of this analytic theory is taken 
from Hubeny, Burrows, & Sudarsky (2003), to which the 
reader is referred for further details. The paper by Cheval- 
lier, Pelkowski, & Rutily (2007) is also of some consider- 
able utility. 

The equation of hydrostatic equilibrium is fundamen- 
tal in atmospheric theory, and obtains as long as the Mach 
number of the gas is low. It can be written: 
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where k is the wavelength-dependent opacity, p p h is the 
mass density at the photosphere, AR c h is the excess radius 
over and above the r p h = § radius (the radius of the tradi- 
tional photosphere), and H is the atmospheric density scale 
height. The latter is given approximately by kT/ pgm p , 
where p, is the mean molecular weight, g is the surface grav- 
ity, T is some representative atmospheric temperature, and 
m p is the proton mass. By definition, and assuming an ex- 
ponential atmosphere, r p h = np p hH = |. For r c h or d to equal 



dP 

- Z = ~ gP 
dP_ _ 
dm ® ' 



(25) 



where P is the pressure, g is the acceleration due to gravity 
(generally assumed constant), p is the mass density, z is the 
altitude, and m is the areal (column) mass density defined 
by dm = —pdz. Note that for a constant g, the pressure 
and the column mass are directly proportional. 

If the atmosphere is convective, the temperature gradient 
follows an adiabat and is given by: 
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dlnP 

where C p is the specific heat at constant pressure and R 
is the gas constant. Given the temperature and pressure at 
any point in the convective region, together with the specific 
heat (which can be determined using the composition), the 
temperature at any atmospheric pressure can be determined. 
With that knowledge, eq. (l26T i can then be used to determine 
the altitude scale. Hence, the combination of the equation 
of hydrostatic equilibrium and eq. d26l i can be used to de- 
termine the adiabatic lapse rate for a given specific heat, 
pressure, and temperature. 

However, for most of a strongly irradiated atmosphere, 
and in radiative zones in general, radiation carries the en- 
ergy flux. The radiative transfer equation is written as 



dm 
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(27) 



where I Vjl is the specific intensity of radiation as a function 
of frequency, v, angle [i (the cosine of the angle of propa- 
gation with respect to the normal to the surface), and the 
geometrical coordinate, taken here as the column (areal) 
mass in. The monochromatic optical depth is defined as 
dr v = Xu dm. 

S v is the source function, given in LTE by 
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(28) 



Xu Xu 

Here k v is the true absorption coefficient, a v is the scat- 
tering coefficient, B v is the black body function, J v is 
the zeroth angular moment of the specific intensity, and 
Xu — k u + a u is the total absorption coefficient. The zeroth 
moment is equal to the radiation energy density, divided by 
An. All coefficients are per unit mass. 

The first moment of the transfer equation is written 

dH v 



dm - X »( Jv Sv) 
which can be rewritten, using equation 

K u {Ju Bit) ■ 



as 
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H v is the first angular moment of the specific intensity and 
is equal to the radiation flux, divided by An. Integrating 
over frequency one obtains 



— = KjJ - K B B . 

dm 



(31) 



where kj and kb are the absorption and Planck mean opac- 
ities, respectively, defined by 



and 
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These two opacities are often assumed to be equal. How- 
ever, one should distinguish them here because the differ- 
ence between kj and k b turns out to be crucial in the case 
of strongly irradiated atmospheres. 

The second moment of the transfer equation is 



dK v 
dm 



XuHu 



and integrating over frequency one obtains 

d^- XHH > 



where 
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which is referred to as the flux mean opacity. 

Finally, the radiative equilibrium equation is written as 



/>oo 

/ k v (J„ - B v ) dv = . 
Jo 



Using the above mean opacities, this can be rewritten 

k.jJ — kbB = . 



(37) 



(38) 



Substituting d38l into (HP , one obtains another form of the 
radiative equilibrium equation, 



dH 



— = , or H = const = (a /An) T e 4 f 



dm 



(39) 



where a is the Stefan-Boltzmann constant. 

From equation d38l . one has B = (kj/kb)J, which 
yields an expression for T through J using the well-known 
relation B = aT A . To determine J, one uses the solution 
of the second moment of the transfer equation K(th) = 
Hth = (a /An) T e 4 ff th, where tr is the optical depth us- 
ing the flux-mean opacity, and expresses the moment K 
through J by means of the Eddington factor, fx = K/ J. 
Similarly, one expresses the surface flux through the sec- 
ond Eddington factor, fg = iJ(0)/J(0) (see also Hubeny 
1990): 
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(40) 



where W is the dilution factor, (R*/a) 2 . This solution is 
exact within LTE and is the generalization of the classical 
Milne atmosphere solution. 

The usual LTE-gray model consists in assuming all the 
mean opacities to be equal to the Rosseland mean opac- 
ity. If one adopts the Eddington approximation (/# = 
1 /3; //j = 1/ v3), then one obtains a simple expression 



Off 
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(41) 



We will consider the most interesting case, namely 
strong irradiation, defined by WT^ ^> T^ s . In this case, the 
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second term in the brackets is negligible, and one may de- 
fine a penetration depth as the optical depth where the usual 
thermal part (oc T c 4 ff ) and the irradiation part (oc WTf) are 
nearly equal, to whit 

Tpcn = W . (42) 

The behavior of the local temperature in the case of a 
strict gray model is very simple - it is essentially constant, 
T = T ee W^T* for r < r pcn , and follows the usual 
distribution T oc r 1 / 4 T e ff in deep layers, r > T pcn . In the 
general case, one has to retain the ratio of the absorption 
and Planck mean (assuming still that the flux mean opac- 
ity is well approximated by the Rosseland mean). In the 
irradiation-dominated layers (r < r pon ), the temperature is 
given by 

T = ~/W 1/4 T^ (43) 

where 

7 ee (kj/k b ) 1/4 ■ (44) 

7 is approximately 1 for no or weakly irradiated atmo- 
spheres. However, in the case of strong irradiation, 7 may 
differ significantly from unity. Moreover, it may be a strong 
function of temperature, and, to a lesser extent, of density. 
This is easily seen by noting that in optically-thin regions, 
the local mean intensity is essentially equal to twice the ir- 
radiation intensity, since the incoming intensity is equal to 
irradiation intensity, and the outgoing intensity is roughly 
equal to it as well. The reason is that in order to conserve 
the total flux when it is much smaller than the partial flux 
in the inward or outward direction, both fluxes should be 
almost equal, and so too must the individual specific and 
mean intensities. 

The local temperature in the upper layers is given, using 
equation (|43l , as 

T/T = 7 (T) . (45) 

It is clear that if 7 exhibits a strongly non-monotonic behav- 
ior in the vicinity of To, for instance if it has a pronounced 
minimum or maximum there, equation d45l l may have two 
or even more solutions. 

This is the origin of thermal inversions, when they occur, 
in the atmospheres of strongly irradiated EGPs (Hubeny, 
Burrows, & Sudarsky 2003). The essential element is the 
differential absorption in the optical on the one hand (since 
this is where most of the irradiating stellar light is found), 
and in the infrared (since this is where most of the emis- 
sion at the temperatures achieved in the atmospheres of hot 
Jupiters occurs). As the above formalism makes plain, an 
inversion is not possible for gray opacities, whatever the de- 
gree of irradiation. For strongly irradiated EGPs, if there is 
not a strong optical absorber at altitude, then stellar optical 
light is absorbed rather deeply in the atmosphere, near pres- 
sures of ~1 bar. This is near and interior to the correspond- 
ing emission photospheres in the near-IR. The result is a 
more-or-less monotonic decrease of temperature with alti- 
tude and decreasing pressure. However, if there is a strong 



optical absorber at altitude, the large value of 7 will allow 
another solution to eq. ( |45T > for which there is an inversion. 
The resulting higher temperatures in the upper atmosphere 
will result in higher fluxes in the mid-IR, e.g. in the IRAC 
and MIPS bands of Spitzer) where their photospheres re- 
side. We see signatures of such inversions in the measured 
spectra of many transiting EGPs, such as HD 209458b, XO- 
1, TrES-4, and HD 149026b. However, to date we do not 
know what chemical species is absorbing at altitude in the 
optical (see §4.3). 

3. LESSONS FROM JUPITER AND SATURN 

Jupiter and Saturn are the largest planets in our solar sys- 
tem and serve as initial paradigms for the atmospheres of 
EGPs. With the largest exoplanets invariably referred to as 
"super- Jupiters," it is instructive to start with an assessment 
of the properties of "regular Jupiters." Both Jupiter and Sat- 
urn have been studied with a variety of remote-sensing tech- 
niques across a wide spectral range, and these approaches 
have provided sufficient information to determine physical 
and chemical properties and their variation in both time and 
space. In the case of Jupiter, these physical and chemical 
properties have been bolstered by in situ observations made 
by the Galileo mission atmospheric probe which both ex- 
tended and served as a measure of ground truth for remote- 
sensing observations. This subsection will concentrate on 
properties of their atmospheres and the spectra of their up- 
welling fluxes in order to provide analogies with extrasolar 
planets which can boast fewer observational constraints. 

Observations of Jupiter and Saturn have been made in 
some detail from the time of Galileo, with serious ground- 
based observations taking place with the advent of photo- 
graphic film and spectroscopy. In the second half of the 
20th century, NASA spacecraft ventured to both - first with 
the Pioneers 10 (Jupiter) and 11 (Jupiter and Saturn) flyby 
spacecraft. Later observations were made with more instru- 
mentation (including an infrared spectrometer, and wide- 
and narrow-angle cameras) by Voyagers 1 and 2, which 
visited both planets, and provided an abundance of infor- 
mation on the dynamics, structure, and composition of both 
Jupiter and Saturn. The Galileo spacecraft orbited Jupiter 
and determined more about the dynamics of individual re- 
gions on the planet and dropped a direct probe into its at- 
mosphere. Most recently, the Cassini spacecraft obtained 
substantial information on the atmosphere, flying by Jupiter 
in a gravity assist on its way to Saturn, which it began or- 
biting in 2005. Most recently the New Horizons spacecraft 
flew by Jupiter on its way to Pluto, collecting information 
on Jupiter's clouds. 

These spacecraft observations were complemented by 
observations from Earth-orbiting platforms beginning with 
observations of auroral emission from Jupiter by the In- 
ternational Ultraviolet Explorer (IUE). They include ultra- 
violet through near-infrared imaging and spectroscopy by 
the instrument complement on board the Hubble Space 
Telescope (HST), disk-averaged spectroscopy by ESA's In- 
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frared Space Observatory (ISO), and even X-ray observa- 
tions by the Chandra Observatory. Additional spectroscopic 
and occultation observations were made using instruments 
on board NASA's Kuiper Airborne Observatory (KAO). 

The ever-improving competence of ground-based obser- 
vational facilities has substantially increased our knowledge 
of these atmospheres. In particular, increasing spatial and 
spectral discrimination became available with larger pri- 
mary mirrors (minimizing diffraction-limited spatial reso- 
lution) and a variety of active optical systems, which have 
been extremely effective in reducing the blurring due to at- 
mospheric turbulence. These have been accompanied at the 
smaller-telescope end by a host of amateur Jupiter watch- 
ers wielding increasingly sophisticated instrumentation, in- 
cluding multi-filtered CCD cameras. 

3.1 Spectra 

Models for the structure of the atmospheres of Jupiter 
and Saturn can be understood in terms of the various re- 
gions in which energy transfer is by either radiation or con- 
vection, similar to stellar atmospheres. To good accuracy, 
much of the atmospheres of both planets can be approxi- 
mated as self-gravitating fluids in hydrostatic equilibrium. 
At depth, the atmospheres of both Jupiter and Saturn are 
dominated by convective processes, and gases rise and fall 
at rates faster than they can radiate away energy. Thus, adi- 
abatic conditions hold (eq. |26| >. 

It is important to note the "principal players" in radia- 
tive transfer in Jupiter and Saturn. In the visible and near 
infrared, the spectrum is dominated by CH4 (methane) ab- 
sorption, with some additional opacity from the H2-H2 
collision-induced fundamental and H2 quadrupole lines at 
higher resolution Figure [5] depicts the near-infrared spec- 
trum of Jupiter from 1 to 6 microns. For Jupiter, the 5-/im 
spectrum is not dominated by reflected sunlight, as is the 
case at the shorter wavelengths, but by thermal emission 
from depth, as the result of a dearth of gaseous absorption 
in this region. It is this atmospheric "window" which al- 
lows glimpses of Jupiter's composition down to a few bars 
of pressure, where H2O vapor was first detected in Jupiter's 
atmosphere from observations using NASA's Kuiper Air- 
borne Observatory, and where high-resolution spectroscopy 
detected PH3 and trace constituents such as ASH3 (arsine) 
and CO (carbon monoxide) at depth. For colder Saturn, 
more of this spectral window is composed of reflected sun- 
light. 

At longer wavelengths, in the middle and far infrared, 
the spectrum is dominated by thermal emission. Figure [6] 
shows the brightness temperature spectrum of both Jupiter 
and Saturn. Rather than flux, these spectra are plotted in 
brightness temperature, which is defined as the temperature 
that a blackbody source would need to emit at a particular 
wavelength in order to match the observed flux. 

Plotting a spectrum in brightness temperature rather than 
directly in flux is a convenient way to display a spectrum 
over a broad range, such as the ones in Fig. [6] because 



the flux varies over several orders of magnitude. In addi- 
tion, plotting the brightness temperature provides a quick 
way to estimate the depth in the atmosphere from which 
most of the radiation emerges. For example, the source of a 
strong emission feature with a high brightness temperature 
is likely to be emerging from molecular emission very high 
in the hot stratospheres of either planet. More details can 
be found in dePater and Lissauer (2004), Irwin (2003), or 
Goody and Yung (1989). 

We also note that Fig. [6] plots its primary spectral scale 
in cm -1 ("wavenumbers") rather than wavelength. His- 
torically, this is the result of such spectra being derived 
from Fourier-transform spectrometers, where inverse of the 
path length difference of a Michaelson-like interferometer 
is the native spectral unit. A wavenumber (cm -1 ) is equal 
to the frequency of the radiation in Hertz (sec" 1 ) divided 
by 29.97. Because the frequency or wavenumber is pro- 
portional to quantum energy, such a plot enables molecular 
lines to be displayed in a manner proportional to the transi- 
tion energy. The result is often (but not always!) a regular 
series of line transitions (or groups of line transitions known 
as manifolds). Some of these are evident in Fig. |6]in the far 
infrared for NH3 (ammonia) or PH3 (phosphine) or near 10 
/im for NH3 . 

Figure [6] shows that the far infrared, where the bulk of 
Jupiter's and Saturn's flux emerges, is controlled by the 
collision-induced absorption of H2 which varies so slowly 
that it constitutes a virtual continuum opacity source. Two 
broad rotational lines are evident, and the far-infrared lines 
of NH3 in Jupiter and PH3 in Saturn show up as absorption 
features on the long-wavelength "translational" component 
of the H2 absorption. When helium collides with H2, a 
slightly different spectrum is created than for pure H2-H2 
absorption, and the translation component is stronger. Fit- 
ting this shape accurately provides a means to determine the 
ratio of He (helium) to H2 in these planets. This part of the 
spectrum also allows one to determine the ratio of ortho-H2, 
where the spins of the hydrogen atoms are parallel, to para- 
H2, where the spins are antiparallel (parallel but pointed in 
opposite directions), because para-H2 alone is responsible 
for all even rotational transitions, such as S(0) and ortho-H2 
alone is responsible for all odd rotational transitions, such 
as S(l). In a quiescent atmosphere, there is an equilibrium 
between these two species of H2 which is purely a func- 
tion of temperature and how their energy levels are popu- 
lated. In reality, however (particularly in Jupiter) H2 gas 
can be transported vertically much faster than this temper- 
ature re-equilibration can take place, and the para-ortho H2 
ratio varies across the face of the planet. This is useful as an 
indirect means to track upwelling and downwellling winds. 
For example, a para-H2 value which is lower than the value 
expected from local thermal equilibrium, but greater than 
0.25, would indicate rapid upwelling from the deep, warm 
atmosphere where its value is close to the high-temperature 
asymptotic value of 25%. However, at the high temper- 
atures of hot- Jupiter atmospheres, transiting or otherwise, 
para- and ortho-hydrogen are in thermal and statistical equi- 
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librium. 

The spectra of Jupiter and Saturn are also filled with dis- 
crete transitions, such as the rotational lines of NH3 and 
PH3 in the far infrared. Lines arising from a combination of 
vibrational and rotational quantum transitions can be iden- 
tified easily in the mid-infrared. Both planets display emis- 
sion features arising from CH4 (methane) and from higher- 
order hydrocarbons, such as C2H2 (acetylene) and C2H@ 
(ethane) which are by-products of the photolysis of methane 
by ultraviolet radiation in the upper atmosphere. They ap- 
pear in emission rather than absorption because the strato- 
spheres of both planets are substantially warmer than the 
upper troposphere. This is the result of warming by sunlight 
absorbed in the near infrared by CH4 and small atmospheric 
particulates. NH3 and PH3 can also be seen as components 
of the mid-infrared spectrum; their lines are in absorption 
because they are most abundant in the troposphere. In fact, 
ammonia is a condensate - much like water vapor in the 
Earth's atmosphere. It does not appear as strongly in Sat- 
urn's spectrum because Saturn is colder than Jupiter, and 
ammonia condenses out deeper in the atmosphere than in 
Jupiter and is at the limits of detectability above the H2 
collision-induced continuum. PH3 lines are, in fact, also 
detectable in Jupiter's spectrum, but are more difficult to 
discern among the forest of NH3 lines in Fig. [6] 

The existence of absorption and emission features from 
H2 and CH4 in the thermal spectrum is also useful for de- 
termining temperatures. Both molecules are well mixed in 
the atmosphere, so any changes which take place in their 
emission from point to point can be attributed to changes 
of temperature rather than abundance. There are techniques 
for inverting the radiative-transfer equations for the emit- 
ted flux, given observations of thermal emission in regions 
dominated by H2 or CH4 which allow temperatures to be 
determined (some authors say "retrieved") as a function 
of altitude. For the H2 absorption features, the relevant 
range is in the upper troposphere - around ~ 10CM-00 mbar 
(~0. 1-0.4 bars) in atmospheric pressure. For CH4 emis- 
sion, the relevant range is the stratosphere at pressures less 
than ^100 mbar. 

Details about molecular transitions in planetary atmo- 
spheres can be found in Goody and Yung (1989) and Irwin 
(2004). Irwin (2004) also discuss in some detail the tech- 
niques for retrieving atmospheric temperatures in the giant 
planets. 

3.2 Clouds 

The visual appearances of Jupiter and Saturn are strongly 
affected by clouds. Jupiter's visible atmosphere is famously 
heterogeneous, and cloud motions and colors have been 
studied for decades. An extensive review of historical and 
amateur images is given by Rogers (1995). Because light 
reflected from a generally homogeneous upper-atmospheric 
haze layer Voyager and Cassini imaging observations of 
Jupiter's and Saturn's cloud fields to determine wind ve- 
locities confirmed and greatly refined ground-based stud- 



ies which showed variable zonal (east-west) winds that 
are variable with latitude, but generally constant in time. 
Jupiter's wind vectors are both prograde and retrograde with 
respect to the rotation of the deep interior, ranging from a 
maximum near 140 m/sec and a minimum near -50 m/sec. 
Doppler tracking of the Galileo probe suggests that these 
zonal winds persist into the atmosphere at least as deep at 
the 20-bar pressure level. Saturn's winds display a strong 
jet of higher speeds, reaching 500 m/sec near the equa- 
tor. A drop of Saturn's maximum jet speed has been de- 
tected in comparing HST-derived winds with those derived 
by Voyager imaging, but it is not clear whether this is due 
to a real change in the wind speed itself or a change of 
the altitude of the cloud particles being tracked. Merid- 
ional (north-south) winds are much smaller and only no- 
table around large, discrete features such as Jupiter's Great 
Red Spot or smaller vortices. Jupiter's largest and longest- 
lived vortices are all anticylonic (rotating counterclockwise 
in the southern hemisphere, clockwise in the northern hemi- 
sphere). Cyclonic storms are also evident and Galileo ob- 
servations noted that they are often the source locations for 
observable lightning; however, their lifetimes are generally 
measured in weeks rather than months to years. 

At near-infrared wavelengths, the dominance of gaseous 
absorption by CH4 and some H2 allows tracking of higher- 
altitude clouds and hazes. Imaging through wavelengths 
that sample particulates at various levels provides a means 
of shaping a 3-dimensional picture of cloud systems. 
Stratospheric hazes can be seen at high latitudes which are 
thought to be products of charged-neutral interactions in 
Jupiter's polar auroral regions. Images of near-infrared re- 
flectivity in Saturn have an overall banded structure which 
is similar to Jupiter's, but Saturn seldom shows individual 
discrete features at these wavelengths that are large enough 
to be detected from the Earth. 

The 5-/im "spectral windows" for Jupiter and Saturn 
provide insight into their cloud systems at several bars of 
atmospheric pressure. The lack of substantial gaseous ab- 
sorption allows cloud tops to be sensed, and thermal emis- 
sion begins to dominate most of the 4.8-5.3 fim region for 
Jupiter and the 5.0-5.3 /j,m region for Saturn. Analogous 
to the visible region, Jupiter's 5-/im appearance is spec- 
tacularly heterogeneous (Fig. |7), with some regions of the 
atmosphere nicknamed "5-/im hot spots" ostensibly cloud- 
less down to deep (5 bars atmospheric pressure) and warm 
(275 K) cloud-top levels. These are mixed with regions 
of both intermediate and cold clouds across the face of the 
disk. The warmer regions are highly correlated with areas 
of darker color in the visible, and the hot spots are corre- 
lated with gray areas. This implies that the lighter colored 
clouds represent an upper layer of tropospheric clouds with 
particles thick enough to be optically thick at wavelengths 
near 5 /im. Such clouds could be related to an ammonia 
ice condensate ("cirrus") cloud near 1 bar to 500 mbar in 
pressure, but spectral signatures of ammonia ice have only 
been seen in particulates undergoing rapid upwelling or re- 
gions of turbulent flow to the northwest of the Great Red 
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Spot. This requires that ammonia condensate to be (i) rim- 
ing nucleation sites made of other materials or (ii) particles 
covered by a coating of other material. The morphology 
of these clouds can also be used to diagnose upwelling and 
downwelling motions, just as the para-H2 distribution dis- 
cussed earlier. Upwelling gas is wet and capable of forming 
clouds, but - after reaching dryer and colder altitudes - the 
condensate has precipitated out and downwelling regions 
will generally be cloudless because of the dry nature of the 
gas. Thus the relatively cold 5-fim interior of the Great Red 
Spot (near the center of Fig. [7) diagnoses a region of up- 
welling, and its ostensibly warm (actually less cloudy) pe- 
riphery indicates regions of downwelling. 

Saturn's thermal emission at ^5 /im reveals a cloud 
structure in the several-bar pressure region every bit as het- 
erogeneous as Jupiter's, although one might say it is a neg- 
ative version. A series of warm zonal bands are interlaced 
with a cold bands, and discrete regions are dominated not by 
clear, warm areas as in Jupiter's atmosphere but by cloudy, 
cold ones. Cassini VIMS observations of this region are 
similar, but take advantage of the entire spectral window 
by observing only thermal emisison on Saturn's night side. 
The persistence of some cloud features that are detectable 
in the visible can also be found in this spectral region, an 
indicator that the dynamics creating them persists over a 
wide vertical range. For example an irregular band near 
high northern latitudes known as the "polar hexagon" attests 
to the presence of meridional waves which perturb the nor- 
mally uniform zonal flow, and the combination of visible, 
near-infrared, 5-/im, and mid-infrared observations shows 
that this feature persists in both the cloud and thermal fields 
from atmospheric pressures less than ^100 mbar down to 
levels with pressures of several bars. 

3.3 Temperatures 

Temperatures in the atmospheres of Jupiter and Saturn 
can be determined using several different techniques. For 
individual locations, occultations of optical light from stars 
by the atmosphere has been used to determine atmospheric 
density as a function of altitude and then transformed into 
temperatures and pressures using the ideal gas law, together 
with the equation of hydrostatic equilibrium and an estimate 
of the mean molecular weight. A similar approach has been 
successful using occultations by the atmosphere of coherent 
spacecraft signals in the radio. 

One such result for a Voyager spacecraft radio occulta- 
tion is shown in Fig. [8] for Saturn, but the Galielo probe 
at Jupiter provided the most direct measurement of tem- 
peratures available for the giant planets (Fig. [8j, with di- 
rect measurements consistent with a "dry" adiabat (i.e., one 
without perturbations to the lapse rate from exchanges of 
latent heat from condensing molecules) from ^400 mbar 
to ^20 bars (thick, asterisked curve in Fig. [8]). At higher 
levels, before the entry velocity of the probe became sub- 
sonic, the deceleration of the probe was Doppler tracked, 
the density of impeding gas determined, and knowledge of 



the mean molecular weight was used to determine the tem- 
perature profile of Jupiter at that location in the stratosphere 
(regular asterisked curve in Fig. [SJ. This profile showed 
the presence of substantial vertical thermal waves in the 
stratosphere, which should provide an substantial additional 
source of mechanical energy for the upper stratosphere. 

The most successful approach for mapping tempera- 
tures across the disk is to invert observations of thermal 
emission in spectral regions dominated by well-mixed con- 
stituents. As described earlier, spectral regions dominated 
by collision-induced absorption by well-mixed H2 and the 
7-/im band of CH4 (Fig. |6]) are used in both planets to map 
temperatures. Fig. [8] contains temperatures derived using 
this technique by the Cassini Composite Infrared Spectrom- 
eter (CIRS) instrument for Jupiter and Saturn near the equa- 
tor and for latitudes of ±30°. This approach has a vertical 
resolution no better than a scale height, and so it acts as 
a low-pass filter for temperature variability with altitude. 
Nonetheless, it provides a good means to derive a range of 
temperatures over time and space, similar to instrumenta- 
tion in weather satellites orbiting the Earth. The technique 
has been expanded from spacecraft spectroscopy to a se- 
ries of middle-infrared filtered images of Jupiter and Sat- 
urn, which provide similar temperature information over a 
wide horizontal range without as much vertical information 
overlap. This has provided a means for ground-based ob- 
servations to supplement the close-up coverage provided by 
spacecraft, with spatial resolutions that have often been as 
good as those available from spacecraft instruments with 
the deployment of mid-infrared imaging instruments on 8- 
meter or larger-class telescopes which minimize the effects 
of blurring by diffraction. The entire host of temperature 
information verifies that the atmosphere is indeed separa- 
ble into largely convective versus largely radiative regions, 
with the convective regions over a variety of different lat- 
itudes and bands converging at depth in both planets (Fig. 

m. 

3.4 Compositions 

Over 99.9% of the composition of both planets by vol- 
ume consists of hydrogen, helium and methane. A sum- 
mary of our current knowledge of the compositions of 
Jupiter and Saturn in the detectable part of the atmosphere is 
given in Table 1 . Their relative abundances, particularly the 
He/H2 ratio, are important constraints on models of their 
formation, evolution and interiors. For Jupiter, the value 
shown in Table 1 is the one derived independently by the 
Galileo Probe Mass Spectrometer (GPMS) and a purpose- 
built Helium Abundance Detector (HAD). The measured 
value for the He abundance may be influenced by its ability 
to form droplets at high pressure and "rain out" of the deep 
atmosphere toward the center of the planet. Neon (Ne) is 
soluble in the He drops and appears to be severely depleted 
in Jupiter with respect to solar abundances. The abundances 
of Ar (argon), Kr (krypton) and Xe (xenon) determined by 
the GPMS are all enhanced by ^2.6 with respect to solar 
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values. 

CH4 is present throughout the atmosphere because it rep- 
resents the simplest stable form of carbon, originally deliv- 
ered to the planet as methane or other carbon-bearing ices 
in planetestimals. This is also true of NH3 for N, H2O for 
O, and possibly PH3 for P, and H2S for S. Ammonia and 
water both condense out in the detectable parts of Jupiter's 
and Saturn's atmospheres; both are responsible for at least 
part of the atmospheric opacity in the submillimeter through 
microwave part of the spectrum and, with phosphine, in 
the 5-/i spectral window. The measurement of O/H via the 
Galileo GPMS value is considered to be too low because 
the Galileo probe descended in a 5-/im "hot spot," that also 
proved to be a region of unusual dessication. Chemical 
equilibrium models predict that H2S should react with NH3 
to form NH4SH (ammonium hydrosulfate) clouds around 
~1— 3 bars pressure (between the upper NH3 cloud and the 
deeper H2O cloud layers. Molecules in chemical disequi- 
librium at the cold upper tropospheric temperatures such as 
PH3 and GeH4 (germane) are present because upwelling 
from warmer depths takes place more rapidly than molecu- 
lar decomposition. 

Methane is present in the stratosphere and is destroyed 
by solar ultraviolet radiation, with higher-order hydrocar- 
bons resulting from the chemical recombination process. 
The most abundant of these are C2H2 (acetylene), C2H4 
(ethylene), and C2H6 (ethane), which act as the principal 
means for stratospheric radiative cooling. Enhanced abun- 
dances of hydrocarbons have been reported in polar strato- 
spheric "hot spot" regions, the result of additional auroral- 
related chemistry. Reports of H2O line emission attest to 
the presence of H2O in the stratosphere, in abundances sug- 
gesting infall from the exterior, either from ring particles or 
from interplanetary ice. 

The isotopic 13 C/ 12 C ratio, determined by the Galileo 
probe for Jupiter and by the Cassini CIRS experiment for 
Saturn are close to the terrestrial value, suggesting a proto- 
solar value with almost no chemical fractionation of carbon 
isotopes in either atmosphere. The observed D/H ratios, de- 
termined from the same set of experiments, also suggests 
protosolar values. In the case of Saturn, this is inconsis- 
tent with interior models which predict an enhancement of 
deuterium because of the larger mass of its core. 

In addition to the reviews already mentioned, detailed 
summaries on several topics concerning Jovian and Satur- 
nian atmospheres, interiors, icy satellites, rings and magne- 
tospheres are given in Bagenal et al. (2004) for Jupiter and 
Brown et al. (2009) for Saturn. 

4. SUMMARY OF EGP FEATURES AND HIGHLIGHTS 

Using the formalisms sketched in §2, and the current 
knowledge of Jovian planet atmospheres summarized in §3, 
one can establish a theoretical framework in which to study 
EGPs for any age, orbit, mass, and elemental composition. 
This theoretical edifice is also understood using the sim- 
ple models provided in §2.4 and §2.6, with which one can 



understand most of their measured and anticipated proper- 
ties. In this section, we highlight a few of the interesting 
ideas and features that have emerged in the last few years 
concerning EGP evolution, generic planet/star flux ratios, 
and the distinctive character of the family of hot Jupiters. 
This summary is in no way comprehensive, but provides 
the reader with a snapshot of some of the salient facts and 
issues in this rapidly evolving field. 

4.1 Evolution of EGP Atmospheres and Chemistry 

Unlike in a star, EGP atmospheric temperatures are suffi- 
ciently low that chemistry is of overriding importance. The 
atmosphere of a gaseous giant planet is the thin outer layer 
of molecules that controls its absorption and emission spec- 
tra and its cooling rate. Molecular hydrogen (H2) is the 
dominant constituent, followed by helium. An EGP's ef- 
fective temperature (T c ff ) can vary from ^2500 K at birth 
for the more massive EGPs or in a steady state for the most 
severely irradiated to ^50 K for the least massive EGPs in 
wide orbits after evolving for billions of years. This wide 
range translates into a wide variety of atmospheric con- 
stituents that for a given mass and elemental composition 
can evolve significantly. 

Without any significant internal sources of energy, after 
formation an EGP gradually cools and shrinks. Its rate of 
cooling can be altered by stellar irradiation, or when old 
and light by hydrogen/helium phase separation (Fortney & 
Hubbard 2004). Jupiter itself is still cooling and its total 
infrared plus optical luminosity is about twice the power 
absorbed from the Sun. Hence, the temperatures and lumi- 
nosities achieved are not just functions of mass and compo- 
sition, but of mass, age, composition, orbital distance, and 
stellar type. 

Since EGPs evolve through atmospheres of various com- 
positions and temperatures, age is a key parameter in their 
study. A theoretical evolutionary scenario for Jupiter itself 
can serve to exemplify these transformations. The evolu- 
tion of other EGPs with different masses and orbits will be 
different in detail, but for the wide-separation variety not in 
kind. The description below is based in part on the review 
by Burrows (2005), to which the reader is referred. Another 
useful review can be found in Marley et al. (2007). 

At birth, Jupiter had a T e g near 600-1000 K and the ap- 
pearance of a T dwarf (Burgasser et al. 2002). It had no 
ammonia or water clouds and, due to the presence of atomic 
sodium in its hot atmosphere, had a magenta color in the 
optical (Burrows et al. 2001). Due to the formation and 
settling to depth of the refractory silicates that condense in 
the temperature range ^1700-2500 K, its atmosphere was 
depleted of calcium, aluminum, silicon, iron, and magne- 
sium 0. Water vapor (steam) was the major molecule con- 
taining oxygen, gaseous methane was the major reservoir 
of carbon, gaseous ammonia and molecular nitrogen were 
the contexts for nitrogen, and sulfur was found in H2S. At 



2 Note that such silicate clouds should exist at depth in the current Jupiter. 
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lower temperatures (<~ 700K), FeS would be the equilib- 
rium reservoir of sulfur. However, as noted above, earlier 
in the planet's history when refractory species condensed 
out and settled gravitationally, the atmosphere was left de- 
pleted of most metals, including iron. The result of this 
"rainout" was that the chemistry simplified and sulfur was 
in the form of H2S, as is observed in the current Jupiter. 
The rainout phenomenon of condensates in the gravitational 
field of EGPs and solar system giants is a universal feature 
of their atmospheres (Fegley & Lodders 1994; Burrows & 
Sharp 1999). 

As Jupiter cooled, the layer of alkali metals was buried 
below the photosphere to higher pressures, but gaseous H2, 
H 2 0, NH 3 , and CH 4 persisted. At a T cff of ~400 K, water 
condensed in the upper atmosphere and water clouds ap- 
peared. This occurred within its first 100 million years. 
Within less than a billion years, when T e g reached ^160 
K, ammonia clouds emerged on top of the water clouds, 
and this layering persists to this day (see section 3). 

Stellar irradiation retards cloud formation, as does a 
large EGP mass, which keeps the EGP hotter longer. Prox- 
imity to a star also keeps the planet hotter longer, intro- 
ducing a significant dependence of its chemistry upon or- 
bital distance. Around a G2V star such as the Sun, at 5 
Gyr and for an EGP mass of 1.0 Mj , water clouds form 
at 1 .5 AU, and ammonia clouds form beyond 4.5 AU (Bur- 
rows, Sudarsky, & Hubeny 2004). Jupiter's and Saturn's 
current effective temperatures are 124.4 K and 95 K, re- 
spectively. Jupiter's orbital distance and age are 5.2 AU 
and 4.6 Gyr. The orbital distance, mass, and radius of a 
coeval Saturn are 9.5 AU, 0.3 Mj , and 0.85 Rj . How- 
ever, as an EGP of whatever mass cools, its atmospheric 
composition evolves through a similar chemical and con- 
densation sequence. Figure [9] depicts the atmospheric tem- 
perature/pressure (T/P) profile for a sequence of 1-Mj , 5- 
Gyr models as a function of orbital distance (0.2 to 15 AU) 
from a G2V star. As the EGP's orbital distance increases, 
its atmospheric temperature at a given age and pressure de- 
creases. Superposed on the plot are the H2O and NH3 con- 
densation lines. Clearly, a given atmospheric composition 
and temperature can result from many combinations of or- 
bital distance, planet mass, stellar type, and age. This lends 
added complexity to the study of EGPs. 

The atmospheres of hot Jupiters at orbital distances of 
^0.02-0.07 AU from a G, F, or K star are heated and main- 
tained at temperatures of ^ 1000-2000 K, roughly indepen- 
dent of planet mass and composition. The transiting EGPs 
discovered to date are examples of such hot objects. At 
high temperatures, carbon is generally in the form of carbon 
monoxide, not methane. As a result of this and the rain- 
out of most metals, EGP atmospheric compositions are pre- 
dominantly H2, He, H2O, Na, K, and CO. There are, how- 
ever, significant day/night differences in composition, tem- 
perature, and spectrum that distinguish a hot Jupiter from a 
lone and isolated planet or star. For instance, on the night 
side, carbon might be found in methane, whereas on the day 
side it might be in CO. This suggests that non-equilibrium 



chemistry, where the chemical rates and the dynamical mo- 
tions compete in determining the composition, might be at 
play. Exotic general circulation models (GCMs) with cred- 
ible dynamics, radiative transfer, chemistry, and frictional 
effects will soon be necessary to understand the equato- 
rial currents, jet streams, day/night differences, terminator 
chemistry, and global wind dynamics of irradiated EGPs, in 
particular, and of orbiting, rotating EGPs, in general. 

There is another interesting aspect to the orbital-distance 
dependence of EGP properties and behavior. At the dis- 
tance of Jupiter, the infrared photosphere is close to the 
radiative-convective boundary, near 0.5 bars. This means 
that a good fraction of the stellar radiation impinging upon 
Jupiter from the Sun is absorbed directly in the convective 
zone, which thereby redistributes this heat more or less uni- 
formly to all latitudes and longitudes of the planet. The up- 
shot of this, and that the luminosity from the core due to the 
remaining residual heat of formation is comparable to the 
stellar irradiation, is that Jupiter (and Saturn) emit isotropi- 
cally in the infrared, the effects of cloud banding and struc- 
ture not withstanding. However, as a gas giant planet moves 
inward towards its primary star, the pressure level of the 
photosphere and the pressure level of the day-side radiative- 
convective zone separate. When the EGP is at ^0.05 AU, 
and after a Gyr, the photospheres have moved little, but the 
radiative-convective boundary is now near a kilobar. More- 
over, the internal flux is miniscule compared with the mag- 
nitude of the irradiation. The result is that heat is not effi- 
ciently redistributed by internal convective motions, but by 
zonal winds in the radiative atmosphere around the planet 
to the night side. This day side/night side dichotomy is a 
central feature of hot Jupiters and can lead to severe ther- 
mal contrasts as a function of longitude. 

Curiously, a star with solar compostion at the edge of the 
light-hydrogen-burning main sequence (M*~75 Mj ) has 
a T c g of ~1700 K. Therefore, an irradiated EGP, with a 
radius comparable to that of such a star, can be as luminous. 

4.2 Planet/Star Flux ratios of Wide-Separation EGPs 

The planet/star flux ratio versus wavelength is the key 
quantity is the study of orbiting planets. Figure [10] taken 
from Burrows, Sudarsky, & Hubeny (2004), depicts orbit- 
averaged (Sudarsky et al. 2000) such flux ratios from 0.5 
/imto 30 ^mfor a l-Mj/5-Gyr EGP in a circular orbit at 
distances of 0.2 to 15 AU from a G2V star. These mod- 
els are the same as those depicted in Fig. [9] The wa- 
ter absorption troughs are in evidence throughout. For the 
closer EGPs at higher atmospheric temperatures, carbon re- 
sides in CO and methane features are weak. For those hot 
Jupiters, the Na-D line at 0.589 /im and the corresponding 
resonance line of K I at 0.77 /im are important absorbers, 
suppressing flux in the visible bands. Otherwise, the opti- 
cal flux is increased by Rayleigh scattering of stellar light. 
As a increases, methane forms and the methane absorption 
features appear in the optical (most of the waviness seen 
in Fig. [TO] for a > 0.5 AU shortward of 1 //m), at ^3.3 
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/im , and at ^7.8 /im . At the same time, Na and K disap- 
pear from the atmosphere and the fluxes from ~1.5 /imto 
~4 /xm drop. For all models, the mid-infrared fluxes long- 
ward of ^4 /im are due to self-emission, not reflection. As 
Fig. [TO] makes clear, for larger orbital distances a separa- 
tion between a reflection component in the optical and an 
emission component in the mid-infrared appears. This sep- 
aration into components is not so straightforward for the 
closer, more massive, or younger family members. For 
these EGPs, either the large residual heat coming from the 
core or the severe irradiation buoys the fluxes from 1 to 4 
/im . The more massive EGPs, or, for a given mass, the 
younger EGPs, have larger J, H, and K band fluxes. As 
a result, these bands are diagnostic of mass and age. For 
EGPs with large orbital distances, the wavelength range 
from 1.5 //mto 4 /im between the reflection and emission 
components may be the least favorable search space, unless 
the planet is massive or young. 

When water or ammonia clouds form, scattering off 
them enhances the optical fluxes, while absorption by them 
suppresses fluxes at longer wavelengths in, for example, 
the 4—5 /im window. Because water and ammonia clouds 
form in the middle of this distance sequence, the geometric 
albedo (A g ) is not a monotonic function of a. This effect 
is demonstrated in Figure QT| taken from Sudarsky et al. 
(2005), which plots the planet/star flux ratio for a Jovian- 
mass planet orbiting a G2V central star as a function of 
orbital distance at 0.55 /im, 0.75 /im, 1 /im, and 1.25 /im. 
Clearly, the planet/star flux ratio does not follow an inverse- 
square law. Clouds can form, evaporate, or be buried as 
the degree of stellar heating varies with varying distance. 
This behavior is included in the models of Fig. [10] but its 
precise manifestations depend upon unknown cloud parti- 
cle size, composition, and patchiness. As a consequence, 
direct spectral measurements might constrain cloud prop- 
erties. As Fig. Q7J suggests, the planet/star contrast ratio is 
better in the mid- to far-infrared, particularly at wide separa- 
tions. For such separations, the contrast ratio in the optical 
can sink to 10~ 10 . 

For the closest-in EGPs (not shown on Fig. ITOb . such as 
HD 189733b, HD209458b, OGLE-TR56b, 51 Peg b, and 
t Boo b, the contrast ratio in the optical should be between 
10~ 5 and 10~ 6 and is more favorable, though still challeng- 
ing. However, for these hottest EGPs, alkali metal (sodium 
and potassium) absorption is expected to dominate in the 
optical and very near-IR. Sodium has already been detected 
in the transit spectrum of HD 209458b. As a result, the ex- 
pectation is that the optical albedos of hot Jupiters should 
be quite low, perhaps less than ~5%. 

Importantly, the transiting EGPs are so hot and their or- 
bits are so favorably inclined that the planet/star flux ratios 
in the mid-IR, in particular in the Spitzer bands, can range 
between ~10 -3 and ~10 -2 . For these objects, the gap 
between the reflection and emission components is com- 
pletely closed. Moreover, the variation in the summed light 
of planet and star during eclipses and the execution of the 
orbit can be used to derive the planetary spectrum itself. 



4.3 Hot Jupiter Highlights 

The discovery of ^60 (and counting) transiting giant 
planets allows one to address their physical structures by 
providing simultaneous radius and mass meansurements. 
However, the proximity of transiting EGPs to their pri- 
maries boosts the planet/star flux contrast ratios in the near- 
and mid-IR to values accessible by Spitzer during secondary 
eclipse. Such measurements in the IRAC (~3.6, ~4.5, 
-5.8, -8 /im) and MIPS (-24 pm) bands, as well as IRS 
spectral measurements from ~5 to ~ 14 /im, are the first di- 
rect detections of planets outside the solar system. In addi- 
tion, Spitzer has yielded photometric light curves as a func- 
tion of orbital phase, and revealed or constrained brightness 
distributions across planet surfaces. This has been done 
(dramatically so) for HD 189733b at 8 and 24 microns, 
for HD 149026b and GJ 436b at 8 microns, and, though 
non-transiting, for both v And b at 24 microns and HD 
179949b at 8 microns. It will soon be done for many more 
hot Jupiters. Excitingly, Spitzer measurements alone have 
yielded exoplanet compositions, temperatures, and longitu- 
dinal temperature variations that have galvanized the astro- 
nomical and planetary communities. 

Furthermore, using NICMOS, STIS, and ACS on HST, 
precision measurements of transit depths as a function of 
wavelength have enabled astronomers to identify atmo- 
spheric compositions at planet terminators. In this way, 
sodium, water, methane, hydrogen, and carbon monoxide 
have been inferred. Moreover, the MOST micro-satellite 
has obtained a stringent upper limit of ~8% to the optical 
albedo of HD 209458b (Rowe et al. 2008). Hence, in ways 
unanticipated just a few years ago, these space telescopes 
are constraining the degree of heat redistribution from the 
day to the night sides by zonal winds (providing a glimpse 
of global climate), are signaling the presence of thermal in- 
versions, and are revealing exoplanet chemistry. The new 
data from Spitzer, HST, and MOST have collectively in- 
augurated the era of remote sensing of exoplanets and the 
techniques articulated in this chapter were developed to in- 
terpret such data. 

As described above, there are many recent highlights in 
the study of giant exoplanet atmospheres, but we will focus 
here on only a few examples. They are the spectral and pho- 
tometric measurements of HD 189733b by Grillmair et al. 
(2008) and Charbonneau et al. (2008) and the photometric 
measurements of HD 209458b by Knutson et al. (2008) and 
Deming et al. (2005), along with their interpretations (see 
also Barman 2008). These objects exemplify the two ba- 
sic classes of EGP atmospheres which have emerged, those 
without and with significant thermal inversions at altitude, 
and are the two best studied transiting EGPs. 

The top panel of Fig. [T2] shows a comparison of the 
HD 189733b data with three theoretical models for the 
planet/star flux ratio. The best fit model (black) is con- 
sistent with most of the data and assumes solar elemental 
abundances and only modest heat redistribution to the night 
side (P„ = 0.15). One of the salient features is the fact that 
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the IRAC 1 (—3.6 pm) to IRAC 2 (—4.5 pm) ratio is greater 
than one. This has been interpreted to mean that the tem- 
perature profile monotonically decreases outward and that 
there is no significant inversion. The left panel of Fig. [13] 
shows such temperature profiles (but for TrES-1), as well 
as the locations of the effective photospheres in the Spitzer 
bands. An appreciable inversion would reverse the IRAC 
1/IRAC 2 ratio, since the effective photosphere in IRAC 2 is 
further out and at much lower pressures than that for IRAC 
1. In addition, the depression of the IRAC 2 flux is con- 
sistent with the presence of carbon monoxide, which has a 
strong absorption feature there. More important is the com- 
parison between theory and the IRS spectrum of Grillmair 
et al. (2008) between —5 and —14 microns. Not only is the 
slope of the spectrum roughly matched, but there is a — 3-cr 
detection of a feature just longward of —6 microns. This is 
best interpreted as the flux peak due to the opacity window 
between the P and R branches of the v-i vibrational bending 
mode of water vapor (see Fig. O . Water would explain 
both the general slope of the IRS spectrum and this peak 
and together they are now taken to be the best indications 
of the presence of water (steam) in an EGP atmosphere. 

The bottom panel of Fig. [T2]compares theoretical mod- 
els for the planet/star flux ratio of HD 209458b at secondary 
eclipse with the corresponding Spitzer photometric data of 
Knutson et. al. (2008) (IRAC) and Deming et al. (2005) 
(MIPS). The differences with HD 189733b are illuminat- 
ing. Here, the IRAC 1/IRAC 2 ratio is less than one and the 
IRAC 3 (-5.8 ^m) flux is higher than the IRAC 4 (-8 pm) 
flux. These features are best interpreted as signatures of a 
thermal inversion at altitude. There must be an absorber of 
optical (and/or near UV) stellar light to create the thermal 
inversion and a hot outer atmosphere (see right panel of Fig. 
[T3l ). with the result that the planet/star ratio is much higher 
in the IRAC 1, 2, and 3 channels than predictions without 
inversions (e.g., the black curve without an absorber). In 
fact, centered near IRAC 3, the IRAC data look like an 
emission feature, centered in the broad water band from 
—4.5 to —8 /im, as opposed to the absorption trough ex- 
pected for an atmosphere with a negative temperature gra- 
dient. Emission features are expected in atmospheres with 
thermal inversions. 

HD 189733b and HD 209458b together represent the 
two classes of hot-Jupiter atmospheres that are emerging 
to be explained. As Fig. [T2] indicates, we can reproduce 
their basic spectral features, but there remain many anoma- 
lies. The fits at —24 pm are problematic and the detailed 
shapes of the spectra are not well reproduced. Moreover, 
the models depend upon the degree of heat redistribution to 
the night side (modeled as P„; Burrows, Budaj, & Hubeny 
2008), which crucially depends upon the unknown zonal 
flows and global climate. More important, though we know 
that there are inversions in the atmospheres of HD 209458b, 
TrES-4, XO-lb, and (likely) HD 149026b, we don't know 
why. What is the nature of this "extra absorber" in the op- 
tical (§2.6)? With what other planetary properties is the 
presence or absence of an inversion correlated? Metallic - 



ity? Gravity? The substellar flux on the planets TrES-1 
and XO-lb are similar, but only the latter shows an inver- 
sion. In the bottom panel of Fig. [12] we have introduced 
an ad hoc optical absorber at altitude merely to determine 
the opacity (n e ) needed to reproduce the measurements. Its 
physical origin is not addressed. Some have suggested it is 
TiO or VO (Fortney et al. 2008). However, these species 
can easily condense out and should not persist for long in 
the upper atmosphere. The so-called cold-trap, which oper- 
ates in the Earth's upper atmosphere to render it dry, should 
also operate for these metal oxides here (Spiegel, Silverio, 
& Burrows 2009). Hence, the nature of the extra optical 
absorber at altitude is currently unknown. One might well 
ask whether the slight thermal inversion introduced due to 
zonal heat redistribution at depth can be the culprit (Hansen 
2008), but the magnitude of the resulting temperature excur- 
sion is slight (< 100 K) and does not approach the — 1000 K 
needed to explain the data. A recent suggestion by Zahnle et 
al. (2009) involves the photolytic production of allotropes 
of sulfur and HS. However, whether a sulfur allotrope or 
any sulfur compound achieves the necessary abundance at 
the requisite upper pressure levels and what its production 
systematics are with the properties of the planet and star 
(metallicity and/or stellar spectrum) remains to be seen. 
The theoretical task is emerging to be complex and chal- 
lenging. 

5. FUTURE PROSPECTS 

The theoretical study of exoplanet atmospheres in gen- 
eral, and of EGP atmospheres in particular, has been ener- 
gized by continuing discoveries and surprises which show 
no sign of abating. In the short term, the last cold cy- 
cle and the upcoming warm cycle(s) of Spitzer will pro- 
vide more light curves and secondary eclipse measure- 
ments. HST will continue to yield high-sensitivity transit 
data. Ground-based extreme adaptive optics systems that 
will enable high-contrast imaging of exoplanets from under 
the glare of their parent stars will be developed and per- 
fected. In the next decade, JWST will constitute a quantum 
leap in sensitivity and broader wavelength coverage, which 
together may lay to rest many of the outstanding questions 
which have emerged during the early heroic phase of EGP 
atmospheric characterization. These include, but are not 
limited to: 

• How do zonal flows in EGP atmospheres affect their 
spectra and phase curves? 

• What is the magnitude of heat redistribution to the 
planet's night side? What is "P„"? 

• How bright is the night side? 

• Upon what do the longitudinal positions of the hot 
and cold spots of strongly irradiated planets depend? 

• What is the depth of heat redistribution? 
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• What causes the strong thermal inversions in a subset 
of the hot Jupiters? 

• Is there a central role for photolysis and non- 
equilibrium chemistry in EGP atmospheres, partic- 
ularly in hot Jupiter atmospheres? 

• Can the chemistry and irradiation be out of phase? Is 
there chemical and thermal hysteresis? 

• How different are the compositions on a planet's day 
side and night side? 

• How different are the temperatures and chemical 
compositions at the terminators of hot Jupiters during 
ingress and egress? 

• Could there be significant gaps in our understanding 
of the dominant chemical constituents of EGP atmo- 
spheres? 

• What is the origin of dissipation in planetary atmo- 
spheres that stabilizes their zonal flows? 

• What are the signatures of non-zero obliquity and 
asynchronous spin in the measured light curves and 
spectra? 

• What are the unique features of the light curves of 
EGPs in highly elliptical orbits and can tidal heating 
contribute to their thermal emissions? 

• What are the thermal and chemical relaxation time 
constants of EGP atmospheres under the variety of 
conditions in which they are found? 

• At what orbital distances are the day-side and night- 
side IR fluxes the same due to efficient heat redistri- 
bution? 

• What is the role of metallicity and planet mass in 
EGP compositions, light curves, and spectra? 

• What are the spectral and evolutionary effects of at- 
mospheric clouds and hazes? 

• Do the convective cores of hot Jupiters cool differ- 
ently on their night and day sides? 

Our theories are poised to improve significantly, with 
credible multi-dimensional simulations undertaken and 
non-equilibrium chemistry and cloud formation addressed. 
Methods are being designed to understand day-night cou- 
pling, the degree of heat redistribution, and the pattern of 
zonal flow. Moreover, with the discovery of more transiting 
planets such as GJ 436b in the ice giant mass and radius 
regime, the variety of signatures will multiply. Finally, 
more wide-separation EGPs will be separately imaged and 
probed. Given the brilliant nature of its past, whatever the 
manifold uncertainties, this field is on a trajectory of fu- 
ture growth and excitement that is establishing it as a core 
discipline of 2 1 st-century astronomy. 
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Table 1 

Compositions of Jupiter and Saturn in Volume Mixing Ratios 



Gas 


Jupiter 


Saturn 


H 2 


0.862±0.003 


0.882±0.024 


He 


0.136±0.003 


0.113±0.024 


CH 4 


2.4±0.5 x 10~ 3 


4.7±0.2 x 10~ 3 


NH 3 


<7xl(T 4 


< 1 x 10~ 4 


H 2 


< 6.5±2.9 x 10~ 5 


0.2-2 x 10 -7 


CO2 




< 3 x 10~ 10 


CO 


1 x 1(T 9 


(1-3) x 10~ 9 


PH 3 


< 6x 1(T 6 


< 4.5 x 10~ 6 


GeH 4 


~2x 10~ 8 


3.4 x 10~ 7 


H 2 S 


1-8 x 1(T 5 


<io~ 6 


C2H2 


< 3 x 10~ 8 


< 1 x 10~ 7 


C2H6 


< 3 x 10~ 6 


<5xl0~ 6 


HCN 


< 1 x KT 10 


< 10- 9 


CH3D 


~2x 1(T 8 


3.2±0.2x 10~ 7 


HD 


-1.1 x 10~ 5 


0.7-1.7 x 10~ 8 


Ar 


1.0 ± 0.4 x 1(T 5 




Ne 


2.3 ± 0.25 x 10~ 5 
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Fig. 1. — The log (base 10) of the mixing fraction as a function of temperature at a total gas pressure of 1 atmosphere for the seven 
molecules shown with solid curves CH4, (black), CO (red), N2, NH3, H2O, H2S, and PH3, and the two molecules shown with 
dashed curves TiO and VO. At 4000 K, CO and N2 are the most stable species, containing nearly all the carbon and nitrogen, 
respectively. With decreasing temperature, CO reacts with H2, forming CH4, which becomes the dominant carbon-bearing species 
at low temperatures, and N2 reacts with H2, forming NH3, which likewise becomes the dominant nitrogen-bearing species at low 
temperatures. Except above about 3000 K, H2O is fully associated containing nearly all the available oxygen that is not bound in CO. 
Below about 1600 K, its abundance temporarily falls slightly due to the condensation of silicates which reduce the available oxygen; 
however, the mixing fraction of H2O then rises again when CO is converted to CH4, which releases the oxygen tied up in CO. Finally, 
at 273 K H2O drops effectively to zero due to the condensation of ice. With decreasing temperature, both TiO and VO rise as they 
associate, then sharply drop to effectively zero when condensates involving Ti and V form. [Taken from Sharp & Burrows (2007)] 
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Fig. 2. — The log (base 10) of the monochromatic absorption a in cm 2 molecule _1 as a function wavelength A in fim in the infrared 
at a temperature of 1600 K and a pressure of 10 atmospheres for the vibration-rotation transitions of H2O, NH3, and CH4. The 
contribution due to different isotopes is included. For this plot we chose a high enough temperature and pressure to ensure that the lines 
were sufficiently broadened to suppress very rapid and large-amplitude fluctuations in the absorption cross section that can otherwise be 
in evidence over short wavelength intervals. In this way, the main band features (which are nevertheless generic for each species) are be 
more easily seen. At significantly lower pressures the broadening of the lines is much smaller and the absorption can change so rapidly 
in short wavelength intervals that the main features do not show up so clearly. As can be seen here, H2O has a strong absorption feature 
just shortward of 3 fim, and CH4 has a strong peak near of 3.3 fim. In the region of 8 fim to 9 fim all three molecules absorb strongly; 
however, between about 10.5 fim and 11 fim NH3 has absorption which is distinctly higher than that of the other two molecules. When 
the combined opacity is calculated, the individual absorptions must be weighted by the abundances. [Taken from Sharp & Burrows 
(2007)] 
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Fig. 3. — Theoretical optical phase functions of 1-Mj , 5-Gyr EGPs ranging in orbital distance from 0.2 AU to 15 AU from a G2V star. 
Near full phase, the phase functions for our baseline models at larger orbital distances peak most strongly. For the cloud-free EGPs at 
smaller orbital distances (0.2 AU, 0.5 AU, and 1 AU), the phase functions are more rounded near full phase. [Taken from Sudarsky et 
al. (2005)] 
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Fig. 4. — The measured visual phase functions for a selection of Solar System objects. A Lambert scattering phase curve, for which 
radiation is scattered isotropically off the surface regardless of its angle of incidence, is shown for comparison. The phase functions of 
the Moon and Mars peak near full phase (the so-called "opposition effect"). A red bandpass Jupiter phase function, taken from Dyudina 
et al. (2004), is also plotted. [Taken from Sudarsky et al. (2005)] 
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Fig. 5. — Near-infrared spectrum of Jupiter taken in a relatively "cold" but reflective region of the atmosphere by the Galileo Near- 
Infrared Mapping Spectrometer (NIMS) in 1996. All radiances except those near 5 fim are sunlight reflected from clouds. Absorption is 
primarily due to CH4 gas with some contributions from a H2 collision-induced fundamental band near 2.2 fim. From Irwin et al. (1998). 
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Fig. 6. — Middle-infrared brightness temperature spectra of Jupiter and Saturn, derived from Voyager and Cassini infrared spectrom- 
eters. Prominent spectral features are identified. Note that the H2 collision-induced absorption is a smooth continuum in the spectrum, 
with the broad rotational transitions S(0) and S(l) identified. 
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Fig. 7. — 4.85-/im image of Jupiter taken at the NASA Infrared Telescope Facility using the facility NSFCam2 instrument on 2008 Sept 
24. The filter used is 0.24 fim wide and illustrates the primarily thermal radiance emerging from Jupiter's multi-layered cloud systems, 
with more radiance emerging from cloud tops in the deepest atmosphere. Jupiter's Great Red Spot can be seen in the lower center of the 
figure as a region whose periphery is defined by dry, downwelling gas which creates a region relatively clear of clouds. 
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Fig. 8. — Temperature profiles for Jupiter and Saturn. Curves in solid lines are derived from Cassini CIRS (infrared) experiment ob- 
servations, representative of near-equatorial (thick solid) and ±30° latitude (regular solid) curves. Galileo probe Atmospheric Structure 
Instrument (ASI) results are shown by the regular solid curve and asterisks for the accelerometer data in the "inverted" stratosphere 
and by the thick solid curve and asterisks for the direct measurements in the troposphere. A Voyager radio occultation curve is shown 
for Saturn by the dashed curve. Note the warmer temperatures in Jupiter than in Saturn, the "inverted" temperatures in the radiatively 
controlled atmosphere at pressures lower than ~100 mbars, vertical waves measured by the Galileo ASI accelerometer experiment, and 
the convergence of the several temperature profiles at different latitudes with depth. 
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Fig. 9. — Profiles of atmospheric temperature (in Kelvin) versus the logarithm base ten of the pressure (in bars) for a family of irradiated 
1-Mj EGPs around a G2V star as a function of orbital distance. Note that the pressure is decreasing along the ordinate, which thereby 
resembles altitude. The orbits are assumed to be circular, the planets are assumed to have a radius of 1 Rj , the effective temperature of 
the inner boundary flux is set equal to 100 K, and the orbital separations vary from 0.2 AU to 15 AU. The intercepts with the dashed lines 
identified with either {NH3} or {H2O} denote the positions where the corresponding clouds form. [Taken from Burrows et al. (2004)] 
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Fig. 10. — Theoretical planet to star flux ratios versus wavelength (in microns) from 0.5 fimto 30 ^tmfor a 1-A/j EGP with an age of 
5 Gyr orbiting a G2V main sequence star similar to the Sun. This figure portrays ratio spectra as a function of orbital diatance from 
0.2 AU to 15 AU. Zero eccentricity is assumed and the planet spectra have been phase-averaged as described in Sudarsky, Burrows, & 
Hubeny (2003). The associated T/P profiles are given in Fig. [9] Table 1 in Burrows, Sudarsky, & Hubeny (2004) lists the modal radii 
for the particles in the water and ammonia clouds. Note that the planet/star flux ratio is most favorable in the mid-infrared. [Taken from 
Burrows et al. (2004)] 
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Fig. 11. — Planet/star flux ratio as a function of orbital distance at 0.55 (im, 0.75 /jm, 1 /im, and 1.25 pirn assuming a G2V central star. 
In each case, the plotted value corresponds to a planet at greatest elongation with an orbital inclination of 80° . Note that the planet/star 
flux ratios do not follow a simple 1/a 2 law. [Taken from Sudarsky et al. (2005)] 
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Fig. 12. — Top: The planet-star flux ratios at secondary eclipse versus wavelength for three models of the atmosphere of HD 209458b 
with inversions and for one model without an extra upper- atmosphere absorber of any kind (lowest curve). The three models with 
stratospheres have different values of P n (= 0.3, 0.4, and 0.5), but are otherwise the same. The old, default model has a P„ of 0.3. 
This figure demonstrates that models with an extra upper-atmosphere absorber in the optical and with P„>0.35 fit the data; the old 
model does not fit. Superposed are the four IRAC points and the MIPS data at 24 /im from Deming et al. (2005) (square block on 
right). Also included, with a question mark beside it, is a tentative update to this 24 /im flux point, provided by Drake Deming (private 
communication). Bottom: Comparison of spectral observations^ith broadband photometry and theoretical models of the dayside 
atmosphere of HD 189733b. The black points show the mean (unweighted) flux ratios spectra from ~5 to ~14 /im, taken by IRS on 
Spitzer and published in Grillmair et al. (2008). The plotted uncertainties reflect the standard error in the mean in each wavelength bin. 
The filled red circles show broadband measurements at 3.6, 4.5, 5.8, 8.0, 16, and 24 fim from Charbonneau et al. (2008). Shown in color 
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Fig. 13. — Dayside (set of higher curves around 10 bars) and nightside (lower set of curves) temperature-pressure profiles for TrES-1 
(left), as a representative of a transiting EGP without a thermal inversion at altitude and for HD 209458b (right), as a representative of a 
transiting EGP with such a thermal inversion. These profiles incorporate the external substellar irradiation/flux and an internal flux for 
the planet corresponding to the temperature of 75 K at are included for HD 209458b. See text in for a discussion. Figures taken from 
Burrows, Budaj, and Hubeny (2008), to which the reader is referred for further details. 
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